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ABSTRACT 

Two-dimensional  equations  of  motion  of  doubly- rota ted  quartz  plates 
for  the  thickness-shear,  flexure,  and  extensional  vibrations  under  in¬ 
plane  initial  stresses  are  employed  to  predict  changes  In  the  fundamental 
thickness-shear  frequencies  due  to  initial  stresses.  Two  types  of  initial 
stresses  are  considered:  (l)  stresses  due  to  a  pair  of  diametral  forces, 
and  (2)  stresses  due  to  steady  accelerations  for  a  three-point  "T"  shaped 
mount  and  a  four-point  "+"  shaped  mount  configurations. 

Force  sensitivity  and  acceleration  sensitivity  coefficients  are 
computed  and  compared  with  experimental  data  and  existing  computed  results. 
For  both  "T"  shaped  and  "+"  shaped  mount  configurations,  mount  orientations 
corresponding  to  the  maximum  and  minimum  of  acceleration  sensitivity  are 
predicted. 


_ INTRODUCTION 


Changes  in  the  resonance  frequencies  of  doubly-rotated,  circular 
quartz  plates  are  investigated  when  plates  are  subject  to  in-plane  initial 
stresses.  Two  types  of  initial  stresses  are  considered:  (1)  stresses  due 
to  a  pair  of  static  diametral  forces,  and  (2)  stresses  due  to  steady 
accelerations  when  the  plate  is  supported  by  a  number  of  ribbon  supports 
attached  to  the  edge  of  the  plate. 

In  a  previous  paper,  a  system  of  six  two-dimensional  equations  of  motion 

of  crystal  plates  was  derived  for  the  vibrations  of  the  flexure,  extension, 

face-shear,  thickness-shear ,  thickness-twist,  and  thickness-stretch  modes 

subject  to  initial  stresses.^  By  studying  the  free  vibrations  of  these 

six  coupled  modes  (without  initial  stresses)  in  a  number  of  doubly-rotated 

cuts  of  quartz  plate  (yxw£  0/<f> ,  with  6  *  33.9°,  0  <  cj)  <  30°),  we  found 

that  in  the  vicinity  of  the  fundamental  thickness-shear  resonances,  only 

2 

thickness-shear,  flexure  and  extensional  modes  are  predominant.  Hence, 
only  these  three  modes  are  retained  in  the  present  study,  and  their  govern¬ 
ing  equations  are  given  in  Section  II.  In  Section  III,  solutions  of  these 
equations  are  derived  by  the  Rayleigh-Schrodinger  method  of  perturbation. 

It  is  seen  that  the  frequency  changes,  in  turn,  depend  on  the  solutions 
of  two  related  problems,  i.e.,  solutions  of  thickness-shear,  flexure,  and 
extensional  vibrations  without  initial  stresses,  and  solutions  of  initial 
stress  problems  due  to  forces  or  accelerations.  These  solutions  are 
given,  respectively,  in  Sections  IV  and  V.  Finally,  the  changes  of  the 
fundamental  thickness-shear  resonances  are  computed  in  Section  VI. 

The  present  study  is  an  extension  of  our  previous  investigations  in  two 
aspects:  (1)  the  stress-strain  relations  when  referred  to  the  plate  axes 


can  now  accommodate  any  doubly- rota ted  cut  of  plate,  and  (2)  the  coupled 
equations  of  motion  for  plate  vibrations  have  been  extended  to  include  the 
extensional  mode  in  addition  to  the  thickness-shear  and  flexural  modes. 

We  note  that  the  solutions  of  initial  stresses  due  to  diametral  forces  and 
due  to  accelerations  are  unified  as  one  by  realizing  that  the  former  is  a 
particular  case  of  the  latter  as  discussed  in  Section  V. 

For  plates  subject  to  diametral  forces,  the  changes  of  resonance 
frequencies,  represented  by  the  force  sensitivity  coefficient  K^,  are 
predicted  as  a  function  of  the  force  orientation  for  various  cut  orientations 
such  as  AT,  (yxwfc) 10°/33.9°»  FC,  IT,  SC,  and  rotated  X-cut.  The  predicted 

A 

results  are  compared  to  the  experimental  values  of  Ballato,  Ballato  and 
Lukaszek,^  and  calculated  values  of  EerNisse"*  by  a  variational  method. 

For  plates  subject  to  accelerations,  a  three-point  "T"-shaped  mount 
and  a  four-point  "+"-shaped  mount  configuration  are  considered.  The  changes 
of  resonance  frequencies,  represented  by  the  acceleration  sensitivity 
coefficient  K  ,  for  an  SC-cut  plate,  are  predicted  as  a  function  of  the 

Q 

direction  of  acceleration  for  various  mounting  or ientations.  Also  values 
of  |K#|  p  *s  proportional  to  the  n2g"  tip  over  results,  are 

computed  as  a  function  of  mounting  orientation  a. 

It  is  found  that  for  the  MT"-shaped  mount  the  minimum  of  acceleration 
sensitivity  occurs  at  orientation  angle  a  *  -15°  and  a  »  75°#  while  the 
maximum  acceleration  sensitivity  occurs  at  a  *  +15°  and  a  ■  *75°.  For 

the  "+"-shaped  mount,  the  minimum  of  |K  |  occurs  at  a  *  0°  and  a  *  A5°» 

a 

while  Its  maximum  occurs  at  a  ■  15®  and  o  *  75° • 
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f  II.  EQUATIONS  OF  MOTION  FOR  CRYSTAL  PLATES  UNDER  INITIAL  STRESSES 

Let  X  ,  Y  ,  and  Z  be  the  crystal lographic  axes  of  a-quartz, 
c  c  c 

among  which  Z  is  the  axis  of  threefold  symmetry  (optical  axis)  and  X 
c  c 

$  is  one  of  the  three  axes  of  twofold  symmetry  (electrical  axis).  The  relations 

between  the  crystallographic  axes  and  the  plate  axes  of  a  rotated  Y-cut 
and  a  doubly- rota ted  cut  of  quartz  plate  are  shown  in  Fig.  I,  in  which  X, 
i.  Y,  Z  or  Xj ,  Xj,  Xj  denote  the  plate  axes  in  a  right-hi  id  coordinate 

system. 

A  system  of  six  equations  of  motion  for  the  incrementa  itions  of 

C  the  six  lowest  modes  under  static  initial  stresses  was  derived  in  a  previous 

paper. '  By  studying  the  free  vibrations  of  these  six  coupled  modes,  i.e., 

extension,  face-shear,  flexure,  thickness-shear,  thickness-twist  and 

•  thickness-stretch,  in  a  doubly-rotated  quartz  strip  with  a  pair  of  traction- 

free  edges,  we  found  that  for  a  series  of  doubly- rota ted  cuts  with  6  *  33.9° 

but  <J>  -  0,  10°,  15°,  19.1°,  21.9°,  30°,  and  in  the  vicinity  of  thlckness- 

ft  shear  cut-off  frequencies,  the  amplitudes  of  the  thickness-shear,  flexure, 

2 

and  extension  are  much  more  predominant  than  the  other  three  modes. 

Hence  by  neglecting  the  coupling  with  face-shear,  thickness-twist  and 
ft  thickness-stretch  modes  in  Eqs.  (55)  of  Ref.  I,  we  have  the  two-dimensional 

stress  equations  of  motion  for  thickness-shear,  flexure  and  extensional 
vibrations  under  initial  stresses: 

< 


r' 


- 


2  k3ft--(l) 

y  b  pUj 


h;  [w  ♦  ^yso)  ♦  (•  ♦  “i,>)*r)  ♦4:m,)  * 

*  u<:>t<-> . .  T co» 0<o)  t  T(o.u<o,j .  2bpa,o,  , 


f-  [('  *  ^°il«,(0)  *  u{0>t<0>  4  T'»>u'°>  4  T<0'U;°>1 

3xj  M  1 , 1 J  1  1 ,3  5  1  l  ,1  5  1,3  1 

+  !-[[,+  u(0)]t(0)  +  y(0)  (0)  +  T(o,,.<o)  +  (o)  (0)1 

9x3  ('  +  ul,lJt5  +Wl,3t3  5  1,1  3  1,3J 


where  T  and  uj°\  are,  respectively, 

the  non-vanishing  components  of  stress  and  displacement  due  to  the  in- 
plane  initial  stresses,  t^  and  t^  are  the  zero  and  first  order 
incremental  stress  components,  and  u^  and  are  the 

components  of  displacement,  respectively,  associated  to  the  thickness- 
shear,  flexure  and  extensional  modes  of  incremental  vibrations. 

Also  from  Ref.  1,  we  have  the  linear  stress-strain  and  straln- 
dispiacement  relations  for  initial  fields: 


T(0)  - 


2bC  V 


m 


[r 


[%] 


M 


j] 


In  (2).  b  is  one-half  of  the  thickness  and  C  are  the  second  order 

pq 

elastic  stiffnesses  of  the  plate.  In  obtaining  the  last  relation  in  (3)# 
we  have  assumed  that  the  in-plane  rotation  associated  with  initial  deforma¬ 
tion  is  negl igi ble. 

The  stress-strain  and  strain-displacement  relations  for  the  incremental 
vibrations  are  deduced,  respectively,  from  Eqs.  (50)  and  Eqs.  (53)  of 
Ref,  1  as  follows. 

t(0)  -  2bc(0).,  ,n(0) 

p  pq  (p)  (q)  q 

(h) 

c<0>n(,)  . 

P  3  pq  q 


,  ♦  u(0)]u(0)  ♦  U(0)u(0) 

U1 . 1 J  ul . 1  u2,lu2,l 


U(0)u{,) 

u1,3  l 

U(0)u(0)  ♦  U(0)u(0) 

Ul,3 1,1  U2,3U2.) 


•  *  *  ('  ♦ 


.  .  fi  *  u<°nu!'> 


‘tin? 
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where 


-  C  +  C 

pq  pq  pqr  r 


(6) 


and  C  are  the  third-order  elastic  stiffnesses  referred  to  the  plate 
pqr 


axes  (x.)  of  the  crystal  plate.  In  (4),  correction  factors  </  \  have 


r  -------  '  •  (p) 

been  introduced  in  order  to  improve  the  accuracy  of  the  dispersion  relations 


from  the  present  approximate  equations  of  motion.  Their  values  will  be 
given  in  Section  IV  when  the  free  vibrations  of  the  crystal  plates  without 
initial  stresses  are  studied. 

By  considering  incremental  vibrational  modes  dependent  on  Xj  only, 
we  may  wr i te 


(1)  (  \  iwt 

Uj  «  w(Xj)e 


u<°>  =  bv(xj )e 


fait 


(7) 


(0)  _  L  /  ^iwt 


u|  =  bu(x^)e 


Substitution  of  (7)  into  (5),  then  into  (k)  and  (1)  leads  to  the  displace¬ 
ment  equations  of  motion  for  steady  vibration 


fiw,n  +  V,i  +  V  +  V,11  +  V,1  +  V.n  +  V.i  "  'n  w 


'Flw,ll  +  F2W,  I  +  F3w  +  V.ll  +  V,l  +  P6U,»I  +  F7U,1  •  V 


Flw,ll  +  F2W,  1  +  F3W  +  E**V,I1  +  F5V,1  +  E6U,  1 1  +  E7U,1  u 


(8) 


/^2J/2 


where  ft  ■  u)/w^  is  a  dimensionless  frequency,  *  (3<  C^/pb  ) 


2  2 

the  lowest  thickness-shear  cut-off  frequency  and  k  *  *rr  / 1 2 . 


In  (8) ,  Tj,  F.t  and  Ej  are  functions  dependent  on  the  initial 
fields  of  stress  and  displacement  and  the  material  properties  C  and 
C^r.  The  explicit  expressions  of  these  functions  are  given  as  follows 

♦<!)♦«] 

}K  l66 

[c<«(i  . 

f2  ■  TJT  (it"  ('  *  *  cfs>uf°3]  ('  *  um)  |.i 
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+  -!!!)]  (--{!!) 


T? 
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* 
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+  IC 


+  K 


0) 
1 
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Lc5>  2,1  55 


*  [Ul| 
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\ 
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11 
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Once  the  initial  fields  of  stress  and  displacement  are  determined  from 
the  solutions  of  initial  stress  problems  which  will  be  discussed  in 
Section  V,  then  these  functions  of  (9)  may  be  regarded  as  known  functions 
of  Xj  and  x^.  Each  of  these  functions  may  be  separated  into  two  parts  as 


Ts  ■  T,  +  t.  • 

i  i  i  * 


Fi  *  V 


Ei  *  V 


where  T.f  F.  and  E .  are  associated  with  vibrational  motions  without 
initial  stresses,  while  t.,  f^,  and  are  the  quantities  contributed 

by  the  initial  fields.^  Therefore,  by  setting  initial  stresses  and  displace 
ment  gradients  in  (9)  equal  to  zero,  we  obtain 


T2  -  T4  -  T6  -  F1  -  -  FS  "  F7  -  El  -  E3  "  ES  "  E7  -  0  '  ("> 

The  explicit  expressions  for  t.,  f.  and  e.  are  obtained  by  inserting 
(9)  and  (11)  into  (10) . 


I » 
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III.  SOLUTIONS  BY  perturbatioh  method 
For  the  thickness-shear.  Flexure,  and  extensiona)  vibrations  under 
static  initial  stresses,  (8)  can  be  expressed  In  the  matrix  form 


LV  -  XV 


(12) 


2 

where  X  •  -Q, 

is  the  eigenvalue 

,  V  «  (w,v,u 

i)^  the  displacement  vector. 

and  L  the  linear  differential  operator.  By 

inserting 

(10) 

into  (8), 

we 

l 

may  separate  L 

into  two  parts 

L  -  L0  +  Q 

(13) 

% 

where 

i 

i 

~V,1  *  T3 

Vi 

Vi  " 

Lo  * 

F29l 

Vn 

Vn 

« 

_E231 

E43.. 

Vn. 

9 

(14) 

A 

t43ll+t5ai 

Vi 

- 

¥ 

Q  - 

Vn*f2Vf3 

Vl  l+f53l 

f6311+f73l 

e23,  *  e3 

e4311+e591 

e6311+e731 

05) 

Lq  is  the  part  of  the  operator  associated  with  vibrations  without 
initial  stresses,  while  Q  is  the  part  of  the  operator  containing  all  the 
effects  of  the  initial  stresses  inside  the  plate.  In  general,  initial 
stresses  also  affects  the  traction  boundary  conditions,  for  instance,  as 
shown  in  Eqs.  (52)  of  Ref.  1.  However,  we  will  neglect  contributions  of 
initial  stresses  to  the  traction  boundary  conditions  by  assuming  that 

< 

these  effects  at  the  edge  of  the  plate  have  little  Influence  on  the 


O 
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1  i  ® 


fundamental  thickness- shear  vibration  whose  displacement  is  mostly  confined 
in  the  central  portion  of  the  plate.  For  instance,  by  neglecting  the 
initial  fields  in  (4),  the  traction-free  conditions 


t{'>  .  ,<°>  -  t<°>  .  0 


at  x,  -  +a 


may  be  expressed  by 


PV  -  0 


at  x,  ■  +a  , 


where  the  linear  operation  P  is  given  by 


C.l9l 

kk6C66 

kk6C16 


kk6C6691  k,CIC1631 
kk6C1691  kk1C1I91 


From  previous  calculations  of  the  Initial  stress  fields,  *  we  found 
that  the  values  of  t.,  f.,  and  e.  are  of  several  orders  of  magnitude 

ill 

smaller  than  those  of  T.,  F. ,  and  E.,  respectively.  Hence,  we  employ 
the  Rayleigh-Schrodinger  method  of  perturbation  to  calculate  the  frequency 
changes  due  to  initial  stresses. 

Let  and  A^  be  the  eigenvector  and  eigenvalue,  respectively, 

of  the  nt*1  mode  satisfying  (12)  and  (17).  Therefore, 

IV  -  A  V  ,  PV  *  0  at  x.  »  +a  (19) 

n  n  n  n  ] 


whe  re 


L  -  LQ  +  eQj  ,  Q  -  eQ1  .  (2 

A  small  parameter  e  is  introduced  in  (20)  to  indicate  the  smallness  of 
Q  as  compared  to  Lq,  while  Qj  is  assumed  to  have  the  same  order  of 


magnitude  as  Lq.  We  let 


i 


Ll 


e' 


Substitution  of  (21)  into  (19)  and  collection  of  the  zero  and  first 
order  terms  of  e  leads  to 


L0Vn0  ‘  ^nOVnO’  PVnO  '  0  •<*,-« 


lOVnl  *  ^lVnO  '  Wnl  *  XnlVn0‘  p»„|  *  0 

We  see  that  (22)  corresponds  to  plate  vibrations  without  initial 
stresses  and  with  traction-free  edge  conditions.  Lq  is  a  symmetric 
operator  and  V^'s  are  free-vibrational  modes  forming  a  complete, 
othonorma!  set  (after  normalization).  Therefore 


V  .  •  V  -  -  6 
mO  nO  mn 


In  (23),  X  .  and  V  .  are  the  first  order  perturbations  of  X 
n  l  n  I  n 

and  V^,  respectively.  We  let 

V  .  «  I  a  V  A  { 

nl  L  nm  nO 
m 

By  inserting  (25)  into  the  first  equation  of  (23),  multiplying  the 
resulting  equation  by  VnQ  and  employing  (24),  we  obtain 

•  »„o  •  «l  •  Vn0 


or  eX  ,  -  V  _  •  eQ,  •  V  *  »  V  >  •  Q  •  V  . 

nl  nO  I  nO  nO  nO 

The  second  equation  of  (23)  is  also  satisfied  by  (25)  because 
satisfy  the  second  of  (22). 

We  choose  n*t  to  indicate  the  fundamental  thickness-shear  mode 


of  vibration.  From  the  second  of  (21),  we  have 
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IV.  THICKNESS-SHEAR,  FLEXURAL  AND  EXTENSIONAL  VIBRATIONS 

OF  DOUBLY-ROTATED  QUARTZ  PLATES 

2 

For  the  resonance  frequencies  (A^q  »  -J2jq)  an<*  vibrational 
modes  *  (w,  v,  u)\  we  should  look  for  solutions  satisfying  the  two- 

dimensional  governing  plate  equations  and  traction-free  conditions  at  the 
edge  of  a  circular  crystal  plate.  An  analytical  solution  for  such  a 
problem  is  extremely  complex  and  difficult  to  obtain,  if  not  impossible. 
However  at  predominantly  fundamental  thickness-shear  resonances,  the 
vibrational  motion  is  mostly  confined  in  the  middle  portion  of  a  flat 
plate.  This  mode  is  concentrated  even  more  in  the  central  portion  of  a 
contoured  or  electrode  plated  plate. Furthermore,  the  thickness- 
shear  resonance  frequencies  are  shown  to  be  insensitive  to  the  width  of 

g 

a  rectangular  plate  by  Sykes  and  predictions  from  one-dimensional 
solutions  are  shown  to  be  very  accurate  when  they  are  compared  with  expert  - 

g  t  q 

mental  data.  Therefore  in  this  section  of  the  paper,  the  circular 

plate  is  approximated  by  a  square  plate  of  equivalent  area  and  one-dimensional 
solutions  are  obtained  to  predict  resonance  frequencies  and  vibrational 
modes . 

The  one-dimensional  equations  of  motion  and  tract ion-f ree  edge  con¬ 
ditions  at  Xj  ■  +a  (2a  */tT~  R,  R  *  radius  of  the  circular  plate)  are 
given  by  (22),  with  LQ  and  P  defined,  respeci tvely,  by  (14)  and  (18), 
and  n  *  l . 

We  assume  that  displacement  \l^  *  (w,v,u)^  has  the  vibrational  form 
Alq  cos  Vl 
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l  A2q  Sin 

q-i 


l  A,  sin  £  x 
L  3q  n 


q-J 


«  1 
q  1 


(29) 


By  substituting  (29)  into  the  first  of  (22),  we  have  the  dispersion 
relation  (ft  vs.  q*l,2,3) 


“Vl6 


'  ^16 
3fiH2566 
"^66 


•  Vl6 
'  V66.2. 

n2_c  .  5all_L 

W  C66  ,.2 


3< 


(30) 


where  E  *  £  .  are  the  dimensionless  wave  numbers  and  the  roots  of  the 

q  qb 

above  bicubic  equation  in  ft  and  \  for  a  given  value  of  ft,  and 


C  *  ic /  /k  C// 

pq  (p)  (q)  pq  66 


(3D 


From  the  linear,  homogeneous  equations  governing  the  displacement 
amplitudes,  we  obtain  the  amplitude  ratios 


A.  :  A.  :  A,  *  1  :  ou  :  a, 
lq  2q  3q  2q  3q 


q-1.2,3 


(32) 


where 


a,  *  D 
2q 


SC66 


-«?,6 

3^q566 


a,  ■  D 

3q 


-1 


"^16 


^qC66 
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Correction  factors  introduced  in  (4)  and  consequently  appeared 

in  (9).  (11),  (30)  and  (31)  are  given  below 

1  if  p  -  3,5 

X  “  —— —  if  p  «  2 ,k 

/T2 

X,  .  »  4  (3 It) 

(p)  K,  if  P  *  1 


-  ft  i  f  p  *  6 


XK,  •  -  ft 

6  m  c 


In  a  previous  study  of  free  vibrations  of  the  first  six  modes  in 
doubly-rotated  quartz  plates,  correction  factors  are  introduced  to  correct 
the  cut-off  frequencies  of  th i ckness-s t retch,  th i ckness- twi st ,  and  thickness 
shear  modes  (or  called  a,  b,  and  c  modes,  respectively)  by  letting 

1  if  p  *  1,3,5 

P  k  «  — —  if  p  *  2,^,6. 


Since  only  three  of  the  six  modes  are  included  in  the  present  study,  two 
new  correction  factors  and  k ^  are  introduced  in  (3*0  to  correct 

the  behavior  of  the  dispersion  relation  from  the  three-mode  theory  as 
compared  to  that  of  the  six-mode  theory,  is  employed  for  correcting 

the  thickness-shear  cut-off  frequencies  and  Kj  for  adjusting  the  slope 
of  the  extensional  branch  (Ej).  The  values  of  k ^  and  for  various 

cuts  of  quartz  plate  are  listed  in  Table  ). 


Dispersion  curves  computed  from  (30)  for  an  SC-cut  quartz  plate  are 

shown  in  Fig.  2  by  the  solid  lines  as  compared  with  those  in  dashed  lines 

2 

computed  from  the  six  coupled  equations  of  motion.  It  may  be  seen  that 
the  three  corresponding  frequency  branches  from  the  three-mode  and  six¬ 
mode  theories  agree  closely  for  frequencies  jp  to  and  including  the 
fundamental  thickness-shear  cut-off  frequency. 

Substitution  of  (29)  into  the  second  of  (22)  for  traction-free  edge 

conditions  leads  to  three  linear  homogeneous  equations  on  A. 

Iq 

fMpq]{Alq}  "  °*  P,q»  »,2,3  (35) 

where 

Mlq  “  C11*q  sin  V 

M2q  =  tKlC6la3qS  +  K<6C66°  +  a2qV]cos  V 

M3q  =  ^1Clla3q^q  +  KK6C16^  +  a2q^q^COS  ^ qa  .  ^ 

For  nontrivial  solutions  of  (35),  we  have 

det  M  «  0  (37) 

PQ 

which  is  the  frequency  equation  for  the  thi ckness-shear ,  flexural,  and 
extensional  vibrations. 

Calculations  of  the  resonance  frequencies  of  SC-cut  plates  are 

performed  and  illustrated  in  Fig.  3.  By  comparing  Fig.  3  with  Figs. 

and  5,  which  are  the  corresponding  frequency  spectra  calculated  from  the 
2 

six-mode  theory,  we  see  that  the  present  three-mode  theory  gives  very 
good  predictions  for  frequencies  of  thickness-shear,  flexural,  and 


extensional  modes. 


From  (35),  we  can  also  calculate  the  amplitude  ratios 
A11  :  A12  :  Al3  "  }  :  62  :  &3 

if  we  choose  as  the  reference.  Then  the  displacement  field  (29) 

can  be  expressed  in  terms  of  a  and  6  (with  6,  *  l) 

oa  a  I 


I  e  cos  e  x, 

q*l 


A11  i.Vzq  $in  Vl 


y  g  a,  s»n  5  x. 
q-1  «  3<»  «  1 


The  remaining  constant  can  be  determined  from  the  normalization 

condi tion  (2k) ,  i ,e. , 


Vjo  *  V,q  *  (w  v  u}  v  dA  =  l 
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V.  INITIAL  STRESSES  IN  CIRCULAR  PLATES 

Two  types  of  initial  stress  problems  are  considered  in  this  section. 

1 .  Circular  plates  subject  to  a  pair  of  diametral  forces. 

A  circular  plate  of  diameter  d(«  2R)  and  thickness  2b  is  referred 
to  a  rectangular  coordinate  system  X#  Y,  Z  (or  Xj)  with  the  XZ  plane 
being  the  middle  plane  of  the  plate.  Let  N  be  the  magnitude  of  the 
compressive  diametral  forces  and  ip  (the  azimuth  angle)  to  denote  the 
orientation  of  the  force  with  respect  to  the  X  axis  as  shown  in  Fig.  6. 

2.  Circular  plates  subject  to  in-plane,  steady  acceleration. 

The  same  plate  is  now  supported  by  a  number  of  metal  ribbons  attached 
to  the  edge  of  the  plate  as  shown  in  Fig.  7*  The  locations  of  the  supports 
are  denoted  by  angles  a.,  i  *  1,2,3*  **s  for  s  supports  (s  >_  2) .  The 
orientation  of  the  body  force  G  (force  per  unit  volume)  which  has  the 
opposite  sense  to  the  acceleration  is  denoted  by  the  angle  with 
respect  to  the  *  axis  as  shown  in  both  Figs.  1  and  7.  Let  N.  and  T. 
be,  respectively,  the  normal  and  tangential  components  of  the  force  from 
the  support  at  a.  to  the  plate. 

Effects  of  initial  stresses  of  both  types  on  the  changes  in  the 
thickness-shear  resonance  frequencies  had  been  studied  previously  for 
rotated  Y-cuts  of  quartz  plate .  ^  In  both  cases  the  initial  stresses 
were  obtained  from  solutions  for  isotropic  plates,  then  strains  were  cal¬ 
culated  from  the  anisotropic  stress-strain  relations.  The  predicted 
results  were  compared  with  experimental  data  with  good  agreement  for 
AT-  and  Y-cuts. 

In  the  present  paper,  the  same  approach  is  adopted.  Since  the  initial 
stresses  due  to  acceleration,  problem  (2),  have  been  given  in  detail  in  Ref. 

3.  Eqs.  (1)  -  (11),  they  will  not  be  duplicated  here. 


We  would  like  to  note  that  the  solution  for  the  problem  (l)t  Initial 
stresses  due  to  the  diametral  forces,  can  be  deduced  from  that  of  the 
problem  (2)  by  setting 

s  -  2,  G  -  0, 

c»j  -  ty,  a2  -  Oj  +  180°, 

N,  -  N2  -  -N, 

T,  -  T2  - 

in  Eqs.  (1)  -  (11)  of  Ref.  3* 

In  a  study  by  Janiand,  Nissim,  and  Gagnepain,^  initial  stresses 
were  obtained  directly  for  the  anisotropic  plates  and  employed  to  predict 
the  changes  in  frequencies.  It  showed  some  improvement  in  predictions 
of  frequency  changes  due  to  diametral  forces,  but  almost  no  improvement 


for  the  case  of  acceleration  effect. 


VI.  CHANGES  IN  RESONANCE  FREQUENCIES 


2 

Once  the  resonance  frequencies  Xjq(»  ”^10^  's  °*>ta'nec*  ^rom  (37), 
displacement  field  Vjq  from  (39),  and  Q  from  (15)  and  initial  fields 
given  in  Section  V,  the  frequency  changes  due  to  initial  stresses  can  be 
obtained  by  numerical  integration  of  (28) 


Af 


2X 


10 


10 


viodA 


where  A  is  the  electroded  area  of  the  plate. 


(28)' 


(l)  Effect  due  to  a  pair  of  Diametral  Forces 

1 2 

Following  Ratajski's  definition  of  force  sensitivity  coefficient 


we  have 


Af 

f„ 


f07m  -q  nr  v66 


(42) 


where  N  is  the  magnitude  of  diametral  forces,  m  the  order  of  harmonic 
overtones  of  the  thickness-shear  frequencies,  and  0,<t>  are  the  angles  of 
rotation  for  doubly-rotated  cut  as  shown  in  Fig.  1. 

Computations  of  K^.  as  a  function  of  the  orientation  of  the 
diametral  forces  are  made  for  circular  quartz  plates  of  AT-cut  (yx£)33-9 
(yxw£) 10°/33.9°,  FC-cut  (yxw£) l5°/33.9°,  IT-cut  (yxwH) 19°/33.9°, 

SC-cut  (yxw£)220/33«9o*  and  rotated  X-cut  (yxw£) 30°/33.9° .  The  results 


are  shown  in  Figs.  8  to  13*  respectively.  They  are  compared  with  the 

k  5 

experimental  results  of  Ballato,  Ballato  and  lukaszek,  and  the  calculate 
result  of  EerNisse.5  The  results  of  Refs.  5  and  11  both  show  the  Improve¬ 


ment  of  predicted  frequencies  by  employing  the  initial  stress  field  of 
anisotropic  plates  under  a  pair  of  diametral  forces. 
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(2)  Effect  due  to  Steady  Accelerations 

We  define  the  acceleration  sensitivity  coefficient 


K 


a 


Af  8/tt 

^0  / 
ngd/pC66 


F(6, ♦;(»,;*) 


(*3) 


2 

where  B  ■  ngp(7rbd  /2)  Is  the  magnitude  of  body  force  on  the  plate  and 

ng  the  acceleration  in  terms  of  the  acceleration  of  gravity. 

As  indicated  in  (43),  K  is  also  dependent  on  the  positions  of 

a 

supports  of  the  plate.  Two  support  configurations  are  considered  as 
f o l I ows . 

(a)  Three-point  >tTM  shaped  mount 

A  support  conf igurat ion  of  a  circular  plate  with  three  supports 

of  which  two  supports  are  180°  apart  and  the  third  is  90°  from  the 

first  two  forming  a  “TM  shape  is  called  a  nTM  shaped  mount  as  shown  in 

Fig.  14.  The  angle  a  denotes  the  orientation  of  the  "Tn  shaped  mount, 

while  ip  indicates  the  orientation  of  the  body  force  due  to  acceleration. 

Calculations  are  made  for  a  circular  SC-cut  plate  (yxwO 23. 75°/33.9° • 

K  is  computed  as  a  function  of  ip  for  a  *  0°,-15°,-30°,  and  -45°, 
a 

and  is  plotted  in  Fig.  14,  We  see  that  both  the  magnitude  and  location 
of  |K  |  change  as  the  support  orientation  a  changes.  By  repeating 

d  fTl3X 

the  above  process  of  calculation,  we  obtain  I K  |  as  a  function  of 
r  *  1  a 'max 

a  as  given  in  Fig.  15.  It  can  be  seen  from  Fig.  15  that  the  minimum 
acceleration-sensitivity  occurs  at  o  ■  -15°  and  +75°»  while  the  maximum 
accelerat ion-sensi ti vi ty  occurs  at  a  ■  +15°  and  a  ■  -75°*  The  variation 
of  acceleration  sensitivity  as  a  function  of  a  for  "T"  shaped  mount 
has  been  studied  experimentally  by  Goldfrank  and  Warner  and  their  experi¬ 
mental  values  are  shown  as  solid  dots  in  Fig.  15  for  comparison.'^ 
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The  result  shown  In  Fig.  15  Is  also  consistent  qualitatively  with  the 

1  4 

experimental  data  by  Kusters,  Adams  and  Yoshida,  for  their  study  Is 
on  a  two-point  mount  configuration  (see  Fig,  2  of  Ref.  14).  We  note  the 
convention  for  the  mounting  angles  in  Ref.  14  is  the  opposite  of  that 
employed  in  the  present  paper. 

(b)  Four-point  Shaped  Mount 

A  four-point  support  configuration  of  a  circular  plate  for  which 

each  support  is  90°  from  its  neighboring  supports  is  called  a  four-point 

M+"  shaped  mount.  Computations  on  an  SC-cut  circular  plate  have  been 

carried  out  in  similar  manner  as  in  case  (a).  The  plots  of  K  vs.  ^  for 

a 

various  values  of  a  and  |K  I  vs.  a  are  given  in  Figs.  16  and  17, 

'  a 1  max  * 

respectively. 

Figure  17  reveals  that  the  minimum  of  acceleration  sensitivity  occurs 

at  a  *  0°  (or  90°)  and  a  *  45°  (or  -45°)  ,  while  the  maximum  occurs  at 

a  *  15°  and  a  *  75°.  It  may  be  seen  that  there  exists  a  wider  'Valley11 

at  a  *  45°  than  that  at  a  *  0°.  It  therefore  appears  to  suggest  that 

it  is  more  likely  to  obtain  experimentally  a  |K  |  «  minimum  near  a  ■  45° 

a 

than  at  a  *  0°,  Nevertheless,  a  minimum  of  |K  |  at  a  ■  0°  had  been 

3 

attained  exper  imental  ly  by  Warner.^ 

The  effect  of  small  deviations  of  the  support  locations  from  the  90° 
angle  in  the  "+"  shaped  mount  has  also  been  investigated  by  allowing  one  or 
two  supports  having  an  angle  of  deviation  6  for  -3°  6  <  3°*  Results  from 

several  series  of  computations  shows  that  the  magnitude  of  acceleration 
coefficient  |K  |  increases  as  1 6 1  increases  for  small  deviations 
“3°  <_  6  £  3°.  However,  the  |K  |  vs.  a  relation  remains  practically 
the  same  as  that  shown  in  Fig.  17,  therefore  it  is  independent  of  6. 
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In  all  the  calculations,  we  have  employed  Bechmann's  values  of  the 
second-order  elastic  stiffnesses1^  and  those  of  Thurston,  McSkimin  and 
Andreath  for  the  third-order  elastic  stiffnesses1^  of  quartz.  For  calculat¬ 
ing  support  reactions  from  ribbon  supports  due  to  accelerations,  we  have 
assumed  each  ribbon  support  is  made  of  nickel  and  acts  as  a  cantilever 
with  length  i  «  6.35  mm  and  rectangular  cross  section  (h^  *  0.076  mm 

and  h^  ■  1.270  mm).  The  Young's  modulus  for  nickel  is  E  *  4.82  x  10^ 

2 

dyn/cm  .  For  more  detail,  see  Eq.  (2)  of  Ref.  3. 
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Dispersion  curves  of  SC-cut  quartz  plate  for 
waves  propagating  in  the  x.-di rect ion. 


Fig.  9  Kf  vs.  <j/  for  quartz  plate  (yxw®.)  10°/33»9°. 


vs.  ip  for  IT-cut  plate  (yxwd)  19.1V33.90. 


for  SC-cut  plate  (yxw£)  2l.9°/33.9 


Acceleration  sensitivity  coefficient  K  as  a  function  of  the  orientation 
of  body  force  for  SC-cut  plate  with3  ,,TM  shaped  mount. 


vs.  af  for  SC-cut  plate  with  "T"  shaped  mount. 


vs.  Ip,  For  SC-cut  plate  with  "V  shaped  mount. 


vs.  a,  for  SC-cut  plate  with  ,,+"  shaped  mount. 


